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Abstract 

This paper discusses the derivation of an effective shell-model hamiltonian 
starting from a realistic nucleon-nucleon potential by way of perturbation 
theory. More precisely, we present the state of the art of this approach when 
the starting point is the perturbative expansion of the Q-box vertex func- 
tion. Questions arising from diagrammatics, intermediate-states and order- 
by-order convergences, and their dependence on the chosen nucleon-nucleon 
potential, are discussed in detail, and the results of numerical applications for 
the p-shell model space starting from chiral next-to-next-to-next-to-leading 
order potentials are shown. Moreover, an alternative graphical method to de- 
rive the effective hamiltonian, based on the Z-box vertex function recently in- 
troduced by Suzuki et al, is applied to the case of a non-degenerate (0+2) hw 
model space. Finally, our shell-model results are compared with the exact 
ones obtained from no-core shell-model calculations. 
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1. Introduction 



The shell model is the basic theoretical tool for the microscopic descrip- 
tion of nuclear structure. As is well known, this model is based on the 
hypothesis that, as a first approximation, each nucleon inside the nucleus 
moves independently from the others, in a spherically symmetric mean field. 
The nucleons then arrange themselves into groups of levels, the "shells" , well 
separated from each other. 

The shell model reduces the complex nuclear many-body problem to a 
very simplified one, where only a few active nucleons (valence nucleons) in- 
teract in a truncated model space spanned by a single major shell above an 
inert core. 

This scheme does not take into account neither the degrees of freedom of 
the core nucleons nor the excitations of the valence nucleons into the shells 
above the model space. Actually, the physical observables must be calculated 
within an effective theory, namely the shell-model hamiltonian has to take 
into account all the degrees of freedom that are not considered explicitly. 

The derivation of the effective shell-model hamiltonian may follow two 
paths. The first and very successful one is phenomenological, where the 
one- and two-body components of the hamiltonian - the single-particle (SP) 
energies and the residual interaction - contain free parameters that are fixed 
to reproduce a set of experimental data. This can be done either using an 
analytical expression for the residual interaction with adjustable parameters, 
or treating the hamiltonian matrix elements directly as free parameters (see 

Q0). 

The alternative approach to the derivation of the effective shell-model 
hamiltonian is the microscopic one, where one starts from a realistic nucleon- 
nucleon (NN) potential, and possibly a three- nucleon one (3iV), and derive 
the effective hamiltonian in the framework of the many-body theory. This 
means that when diagonalizing the shell-model hamiltonian in the model 
space, its eigenvalues belong to the set of eigenvalues of the full nuclear 
hamiltonian, defined in the whole Hilbert space. 

The most successful way to derive a realistic effective shell-model hamilto- 
nian is rooted in the energy-independent linked-diagram perturbation theory 
jsj], which has been widely employed in shell-model calculations during the 
last forty years (see review papers 0, lEf). The core of this approach is the 
perturbative expansion of a vertex function, the so-called Q-box, as a col- 
lection of irreducible valence-linked Goldstone diagrams. The Q-box is then 
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employed to solve non-linear matrix equations to derive the desired effective 
hamiltonian, which can be done by way of iterative techniques jfjj]. In all ap- 
plications the Lee-Suzuki and Krenciglowa-Kuo iterative methods have been 
commonly employed, both of them based on the calculation of the derivatives 
of the Q-box with respect to the energy 

Recently, Suzuki and coworkers [7| have introduced an alternative way 
to derive the shell-model effective hamiltonian, which is a graphical method 
based on the introduction of another vertex function, the so-called Z-box, 
whose main advantages are that only the first derivative of the Q-box is 
needed, and that it can be easily extended to the case of a non-degenerate 
model space. 

In such a scenario, one should keep in mind that in modern nuclear struc- 
ture calculations it has been evidenced the role of 3N forces, in particular for 
light nuclei with A < 12 js), 0, [To| • Within the framework of the shell- model 
effective hamiltonian theory, the inclusion of 3iV forces yields an effective 
hamiltonian that consists of one- and two-body components, including also 
core-polarization effects arising from the 3N force, and an effective three- 
body interaction that should be explicitly considered in the calculations. 
The derivation of such a hamiltonian, even if not requiring any refinement 
of the theory, is a very hard task, and up to the present a shell-model effec- 
tive hamiltonian derived treating on equal footing both NN and 3N forces 
has never been calculated. In this connection it is worth to cite the work 
by Otsuka et al. [11], where only first-order contributions of the normal- 
ordered two-body parts of 3N forces are taken explicitly into account. On 
these grounds, we confine our study, without any loss of generality, to effec- 
tive shell-model hamiltonians derived starting from purely two-body realistic 
potentials. 

The aim of the present work is to try to assess the state of the art of the 
derivation of the realistic shell-model effective hamiltonian from the pertur- 
bative expansion of the Q-box. In particular, we shall discuss in detail the 
behavior of the perturbative series, with respect to both the dimension of 
the space of the intermediate states and the order-by-order convergence. 

In this context, we will show the results of shell-model calculations for 
p-shell nuclei starting from realistic NN potentials based on the chiral pertur- 



bation theory at next-to-next-to-next-to-leading order (N 3 LO) [l2j, |l3|, Il4 
Moreover, we shall present results of calculations within the (0 + 2) hu psd 
model space, using an effective hamiltonian derived by the graphical Z-box 
method. 
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The paper is organized as follows. In Section 2 we give an outline of 
the derivation of the shell-model effective hamiltonian within a perturba- 
tive approach. Section 3 is devoted to the discussion of the problematics 
concerning the diagrammatics and convergence of the Q-box perturbative 
expansion. This is done showing our results for the p-shell nuclei, using 
both and (0 + 2) fiuj model spaces. In Section 4 we c omp are shell-model 
results obtained starting from the N 3 LO potentia l Il2l _ 13| with the exact 
ones provided by the ab initio no-core shell model (la . lol 
concluding remarks are given in Section 5. 



A summary and 



2. The shell-model effective hamiltonian 

In this section, we describe the formalism of the shell-model effective 
hamiltonian theory. 

Let us consider the Schrodinger equation for the A-nucleon system in the 
whole Hilbert space: 

H\<S! V ) = E V \<S! V ) . (1) 

As mentioned in the Introduction, in the frame of the shell-model ap- 
proach an auxiliary one-body potential U is introduced to break up the nu- 
clear hamiltonian as the sum of an unperturbed one-body term Hq, which 
describes the independent motion of the nucleons, and the interaction hamil- 
tonian Hi. 

We then write H = H + H x with H = J2f=i h o = Y.t=i( T i + u i) and 

//i v,.,o;f v -/',)• 

In the shell model, the nucleus is schematized as an inert core plus n 
interacting valence nucleons moving in the mean field Hq. The SP potential 
U generates an energy spectrum organized in shells. The large energy dif- 
ference between shells enables to define the core as the A — n nucleons that 
fill completely the lowest shells. The SP states accessible to the n valence 
nucleons are then the lowest in energy above the closed core. 

Now, it is possible to define a reduced Hilbert space, the model space, in 
terms of a finite subset of d eigenvectors of Hq 

= [a\al ... al]i\c) , (2) 

where |c) represents the inert core, the subscripts 1, 2, n denote the SP 
valence states and i stands for all the other quantum numbers that specify 
the state. 
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More explicitly, we define the projection operators P and Q = 1 — P, 
which project from the complete Hilbert space onto the model space and 
its complementary space, respectively. They satisfy the properties P 2 = P, 
Q 2 = Q, and PQ = QP = 0. In terms of eigenvectors of Ho, P is defined as 

d 

p = j2mm ■ (3) 

i=l 

Our aim is to reduce the eigenvalue problem of Eq. (JTJ) to the model-space 
eigenvalue problem 

H eS P\* a ) = E a P\* a ) , (4) 

where a = 1, .., d and H e g is defined only in the model space. In other words, 
we are looking for a new hamiltonian H that has the same eigenvalues of 
the A-nucleon system hamiltonian H, and satisfies the decoupling equation 
between the model space P and its complement Q: 

QUP = , (5) 

so that the desired effective hamiltonian is H c g = PUP. 

The new hamiltonian % can be obtained by way of a similarity transfor- 
mation defined in the whole Hilbert space: 

H = X~ 1 HX . (6) 

There is of course an infinite class of transformation operators X that 
satisfy the decoupling equation (jSJ). Suzuki and Lee [<3| considered an operator 
X defined as X = e u . Without loss of generality, they chose u so as to satisfy 
the following properties: 

to = QcoP , (7) 
PcoP = QcoQ = PuQ = . (8) 

Eq. (JTJ) implies that 

to 2 = to 3 = ... = . (9) 

The above equation enables to write X — 1 + u, and consequently it can 
be obtained 
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H cS = PUP = PHP + PHQu . (10) 

The operator u may be obtained by solving the decoupling equation (J5j), 
which may be rewritten in the following form 

QHP + QHQuo - uoPHP - uoPHQuo = . (11) 

The above matrix equation is non-linear and can be easily solved, once 
the hamiltonian H is explicitly known in the whole Hilbert space. However, 
this is not an easy task for nuclei with mass A > 2, and has been recently 



done only for light nuclei within the framework of ab-initio approaches [16 
The standard approach to solve Eq. (|TT|) in a shell-model calculation is 

to introduce a vertex function, the Q-box, that can be then expressed as a 

perturbative expansion. 

Let us consider our model space, which we assume to be degenerate: 

PH P = e P . (12) 

Then, taking into account the decoupling equation §5§, the effective hamil- 
tonian Hf = H c q — PH P can be written in terms of u> 

Hf = PHP - PH P = PHiP + PHiQu . (13) 

We now employ the above identity, the decoupling equation (TTTj) . and 
the properties of H and Hi, to obtain a recursive equation for the effective 
hamiltonian Hf. 

First, since Ho is diagonal, we can write the following identity: 

QHP = QH 1 P + QH P = QH 1 P . (14) 
The decoupling equation (|TT|) can then be rewritten in the following form: 



QHiP+QHQco-uiPHoP+PHiP+PH&u;) = QHiP+QHQu-u(e Q P+Hf) = 

(15) 

Using this form of the decoupling equation, we can write the following 
identity for the operator u: 

-o^^-e^W^* ' (16) 



6 



Finally, inserting Eq. (|T5|) into the identity (ITS]) that defines Hf 1 , we 
obtain a recursive equation: 



Hf(u) = PH 1 P + PH 1 Q- 



1 



-Q#iP - PH X Q- 



1 



e - QHQ ^eo-QHQ 
We define the vertex function Q-box by the following identity: 



Q(e) = PH X P + PH X Q- 



1 



rQPiP 



e - QPQ 

so that the recursive equation (TTTI) can be expressed as 



PfH = Q(e ) - PPiQ- 



(17) 
(18) 

(19) 



e - QHQ 

Suzuki and Lee 0] suggested two possible iterative techniques to solve Eq. 
( fl9l) . both of them based on the calculation of Q-box derivatives. In the next 
subsections, we shall briefly describe these two methods, which have become 
known as the Krenciglowa-Kuo (KK) and the Lee-Suzuki (LS) techniques. 

2.1. The Krenciglowa-Kuo iterative technique 

This iterative approach originates from the observation that Eq. (119j) . 
coupled with the recursive expression for the operator u, Eq. (TTB]) . leads to 
the iterative equation: 



£ 

m=0 



-PH X Q 



e - QHQ 



m+1 



QH X P 



[Hf(u n ^)] m . (20) 



It should be noted that the quantity inside the square brackets in Eq. 
20]) . which from now on we write as Q m (e Q ), is equal to: 



Q m (e ) = -PH X Q 



-1 



e - QHQ 



m+1 



QH X P = — 



ml 



d m Q(r]) 
dr] m 



(21) 



This identity allows to rewrite Eq. (120 j) in the final form: 
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r]=e "1=0 

(22) 

The above approach is known as the Krenciglowa-Kuo iterative method, 
since, if we make the assumption Hl s (co ) = Q(e ), Eq. (122]) can be rewritten 

as 

oo 
i=0 

where 



oo , 

z — ' ml 



m=0 



dt] m 



F = Q(e ) 

Fi = Qi(e )Q(e ) 

F 2 = g 2 (eo)g(eo)g(eo) + gi(eo)gi(eo)g(e ) 

(24) 

This is the well-known folded-diagram expansion of the effective hamil- 



tonian as introduced by Kuo and Krenciglowa in [17|, where the following 
operatorial identity has been demonstrated: 



QiQ = -Q J Q , (25) 

the integral sign representing the so-called folding operation 

2.2. The Lee-Suzuki iterative technique 

Suzuki and Lee [f| suggested another iterative technique, which can be 
introduced by rearranging Eq. f[T§|) so as to have an explicit expression of 
the effective hamiltonian Hf s in terms of the operators uj and Q: 

Hf(u) = (l + PH iQ eo _Q HQ ^ Q(eo) • (26) 
The above equation can be rewritten in an iterative form: 
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tffK)= [1 + PHiQ 



1 



e - QHQ 

and similary for the recursive equation ffT6j) : 

1 „ 1 



W n — Q- 



-QH X P - Q- 



-i 



(27) 



(28) 



e -QHQ^ ^e -QHQ 
The iterative procedure can start from the choice Uq = 0, so we can write: 



Hf( Ul ) = Q(e ) 



Q- 



l 



iQ^P 



e - QHQ 

By doing some algebra, it is possible to demonstrate: 



4M 



-PHtQ- 



1 



i 



e -QHQ^ e -QHQ 

-PHiQ 77777^1 



QH X P 



e - Q#Q 

so that, for the next iteration n = 2, one has: 



(29) 



i 



e - QHQ 



wi Q(e ) 



CJ 2 



l-Qi(eo) 

Q Q^iP - Q 

e - QHQ ^e -QHQ 



u x Hf{u 2 ) . (30) 



Finally, the iterative form of the equation for the effective hamiltonian 
within the LS approach reads: 



n-l 



n-l 



QiM^&nM n 



m=2 



k=n— vn+l 



Q(e ) 



(31) 
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The KK and LS techniques to solve the decoupling equation do not nec- 
essarily provide the same effective hamiltonian. In Ref. jfj it has been 
shown that, when converging, the KK iterative procedure provides an ef- 
fective hamiltonian whose eigenstates have the largest model space overlap, 
while the effective hamiltonian obtained with the LS one has eigenvalues that 
are the lowest in energy. 

It is worth noting that both procedures are constrained to choose an 
unperturbed hamiltonian Hq whose eigenstates belonging to the model space 



are degenerate in energy. In Ref. [19| an alternative approach to the standard 



KK and LS procedures has been proposed, which extends these methods to 
the non-degenerate case by introducing multi-energy Q-boxes. However, this 
approach has proved to be quite complicated for practical applications, the 
only one appeared in the literature being that of Ref. 
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In the following subsection, another method to calculate the shell-model 
effective hamiltonian, recently derived by Suzuki et al. 0], will be outlined. 

2.3. The Z(e) vertex function 

From inspection of Eq. (fT8j) . which defines the Q-box, it can be seen 
that when e approaches one of the eigenvalues of QHQ, Q(e) has some poles, 
which can induce instabilities in the numerical derivation. 

To overcome such these difficulties, Suzuki et al. [7] have recently intro- 
duced an alternative vertex function Z(e) that is defined in terms of Q{t) 
and its first derivative: 



1 



Q(e)-Q 1 (e)(e-e )P 



1-Qi(e 

In Ref. [?} it has been shown that Z(e) satisfies the following equation 



(32) 



e + Z(E a ) P\V a ) = E a P\y a ) 



(a = 1, ..,d) 



(33) 



which means that Hf may be obtained calculating the Z-box for those 
values of the energy, determined self-consistently, that correspond to the 
"true" eigenvalues E a . 

To obtain the E a , let us now consider the eigenvalue problem 



t + Z{t) \<t> k ) =F h (e)\<j> h ) 



(* = l,2,...,d) 



(34) 
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where F k (e) are d eigenvalues that depend on e. The true eigenvalues E a 
may be determined by solving the following d equations 

e = F fc (e), (A; = 1,2, ■■ . (35) 

Before discussing how to solve equations (I34II35I) . let us point out some 
interesting properties of Z(e) and of the associated functions F k (e). 

In the neighborhood of the poles of Q(e), the behavior of Z(e) is dom- 
inated by Qi(e) and so it can be written Z(e) « (e — eq)P. This means 
that Z(e) has no poles and therefore the functions F k (e) are continuous and 
differentiate for any value of e. 

Eqs. (I3"5"|) may have spurious solutions, i.e. solutions that do not corre- 
spond to the true eigenvalues E a . However, in Ref. |7| it has been shown 
that the energy derivative of F k (e) becomes zero at e = E a , so the study of 
this derivative gives a criterion to get rid of the spurious solutions. In order 
to solve Eqs. (134"]) and (1351) . so as to derive the effective interaction, one may 
resort to both iterative and non-iterative methods. 

In the present paper we have employed a graphical non-iterative method 
to solve Eqs. 0351) . which we shall now describe. 

As shown previously, the Ffc(e)'s are continuous functions of the energy, 
therefore one of the well-known algorithms to solve nonlinear equations may 
be employed to determine the solutions of Eqs. fl35l) as the intersections of 
the graphs y = e and y = F k (e). 

More precisely, if we define the functions f k (e) as /^(e) = F k (e) — e, the 
solutions of Eqs. (135]) can obtained by finding the roots of the equations 
/fc(e) = 0. To this end, from inspection of the graphs y = e and y = F k (e), 
we determine for each intersection a small surrounding interval [e a , e&] such 
that fk(^a)fk{^b) < 0. The assumption that f k (e) is monotone in this interval 
implies the existence of a unique root, that we can determine very accurately 



employing the secant method algorithm (see for instance Ref. 211]) . 

Once we have determined the true eigenvalues E a , the effective hamilto- 
nian Hf is built up as 

d 

Hf = J2z(E a )\^ a )(4> a \ , (36) 

a=l 

where \<p a ) is the eigenvector obtained from Eq. while (<f> a \ is the corre- 
spondent biorthogonal state ({<fi a \<j>a') — 8aa')- 
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In conclusion, it is worth pointing out that in the previous discussion 
we have considered the case of a degenerate unperturbed model space (i.e., 
PHqP = €qP). However, the above formalism can be easily generalized to 
the non-degenerate case replacing e^P with PH P in Eqs. fl32|) and fl34l) . 



3. The perturbative approach to the shell-model H eff 

The present section is devoted to discuss details and problems of the 
derivation of a shell-model effective hamiltonian by way of the perturbative 
approach. More precisely, the Q-box is calculated perturbatively and then 
employed to derive the effective hamiltonian within the KK, LS, or Z-box 
graphical approaches. It should be pointed out that the above techniques 
lead to effective hamiltonians whose matrix elements differ at most by a few 
keV when considering the KK and LS methods, while when using the Z- 
box method the differences do not exceed 80 keV. This can be seen from 
inspection of Tables IA.3I and IA.4I that can be found in |Appendix A 



3.1. The diagrammatic expansion of the Q-box 

The methods presented in the previous subsections are based on the cal- 
culation of the Q-box function: 

Q(e) = PH X P + PH l Q 7 -^ w ^ QH l P . (37) 
The term l/(e — QHQ) can be expanded as a power series 

QH t Q 



e - QHQ ^ e - QH Q \e - QH Q 

giving rise to a perturbative expression for the Q-box. The diagrammatic 
representation of this perturbative expansion is a collection of diagrams that 
have at least one ifi-vertex, are irreducible (i.e., with at least one line not 
belonging to the model space between two successive vertices), and valence 
linked (i.e., are linked to at least one external valence line) [221 ] . 

Currently, realistic shell-model effective hamiltonians are derived for sys- 
tems with one and two valence nucleons. The former provides the theoretical 
effective SP energies, while the two-body residual interaction V eS is obtained 
from the H eS of the two-valence-nucleon system using a subtraction proce- 
dure [23[. Up-to-date applications include in the Q-box at most Goldstone 
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diagrams up to third order in Hi, which take into account up to 3p — 2h 
excitations for the one valence-nucleon system, and up to 4p — 2h excita- 
tions for the two valence-nucleon system. A comprehensive work concerning 
the evaluation of the linked Goldstone diagrams in an angular momentum 
coupled representation may be found in Ref. 



24 



J 


j 


j 






h 








j 


j 


j 






(a) 



(b) 



Figure 1: One-body first-order diagram. Graph (o) is the so-called self-energy diagram. 
Graph (6) represents the matrix element of the harmonic oscillator potential U = ^muj 2 r 2 . 



The one-body Q-box diagram at first order in H\ is reported in Fig. [TJ All 



other Q-box diagrams up to third order in H\ can be found in Appendix B 



and are reported in Fig. IB.171 IB. 191 (one-body diagrams), and in Fig. IB. 181 
IB~20l |B~2T1 IB~22l (two-body diagrams). 

From inspection of Fig. [TJ it can be seen that the first-order one-body 
diagram is composed of the so-called self-energy diagram (^-insertion dia- 
gram) minus the auxiliary potential [/-insertion. The [/-insertion diagrams 
arise in the perturbative expansion owing to the presence of the —U term in 
H x . 

The (V — [/)-insertion diagrams turn out to be identically zero when 
employing a self-consistent Hartree-Fock (HF) auxiliary potential 2(|. It is 
worth noting that in most applications the standard choice for the auxil- 
iary potential is the harmonic-oscillator (HO) one, and that the (V — U)- 
insertion diagrams are neglected, assuming that the differences between the 
HO and the HF single-particle wavefunctions are negligible. In Subsection 
13.41 we shall discuss about the contribution of these terms, comparing differ- 
ent effective hamiltonians derived starting from Q-boxes with and without 
(V — [7)-insertion diagrams. 

As we pointed out before, in the existing literature the effective hamilto- 
nians are derived taking into account in the Q-box at most diagrams up to 
the third order, being computationally prohibitive to go to higher-order. 

In order to have a better estimate of the value to which the perturbation 
series should converge, it is helpful to resort to the Pade approximant the- 
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ory [25J, |26], and to calculate the Pade approximant [2|1] of the Q-box, as 



suggested in [27 



[2|1] = V» box + V^ box + V* box (l - , (39) 

where Vq box is the square non-singular matrix representing the nth-order 
contribution to the Q-box in the perturbative expansion. 

In the following subsections we show results which we have obtained cal- 
culating H for two valence nucleons outside 4 He doubly-closed core in the 
p-shell model space, and using two potentials based on the chiral perturba- 
tion theory at next-to-next-to-next-to-leading order. One is the well-known 



N 3 LO potential derived by Entem and Machleidt 12|, |l3|, that is charac 



terized by a smooth gaussian cutoff around 2.5 fm . The other potential, 



dubbed N 3 LOW 14{, has a sharp cutoff with a smaller value A = 2.1 fm 



In both cases, for protons the Coulomb force has been explicitly added. 

3.2. Convergence with respect to the intermediate- state space 

The Q space that enters the definition of the Q-box in Eq. ffl8|) is infinite, 
representing the complement of the model space in the whole Hilbert space. 
This implies that in the calculation of the diagrams composing the Q-box 
one should perform an infinite sum over the intermediate states between 
successive vertices. This is unfeasible, so the space of the intermediate states 
has obviously to be truncated. A well established procedure is to introduce an 
energy truncation, i.e. the intermediate states whose unperturbed excitation 
energy is greater than a fixed value E max are disregarded. E max has to be 
chosen sufficiently large to ensure that the results are almost independent 
from its value. 

In this regard, it is appropriate to mention the papers by Vary et al. 



28j, Kung et al. [29J, and Sommermann et al. [30j, where the convergence 
rate of the sum over intermediate-particle states in the second-order core 
polarization contribution to the effective shell-model interaction was studied, 



using realistic potentials renormalized by way of the Brueckner theory [31 . 

In Figs. [2]and[3j the theoretical energies of the yrast states in 6 Li relative 
to 4 He are reported as a function of E max , expressed in terms of the number 
of oscillator quanta N max . The energies reported in Figs. I2"f3l have been 
calculated using an effective hamiltonian derived from the N 3 LO and N 3 LOW 
potentials respectively, including in the Q-box diagrams up to the second 
order in H\. 
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2 4 6 8 10 12 14 16 18 20 22 

Figure 2: Theoretical energies of 6 Li yrast states relative to 4 He, obtained with the N 3 LO 
potential, as a function of iV max (see text for details). 




From inspection of Fig. [2J it can be seen that convergence is reached with 
the N 3 LO potential when including intermediate states whose unperturbed 
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excitation energy is less than E max = 20 Huj, with ftw = 19 MeV. This value 
of the harmonic oscillator parameter is close to the one provided by the 

- 1 / 3 - 25A- 2 / 3 for A 



expression 
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fkj = 45A~ L ^ - 25A-'^ for A = 4. 
A faster convergence can be obtained using the N 3 LOW potential, since, 
as shown in Fig. |3l the energies are practically stable from E max = 10 hw 
on. This is related to the fact that the two potentials are characterized by 
largely different cutoffs, N 3 LOW being a low-momentum potential. 

In most applications a subtraction procedure [23| is used so as to retain 
from Hf 1 only the effective two-body interaction V , while the SP energies 
are taken from experiment. 

We have therefore found it worthwhile to study the convergence properties 
of V eG when enlarging the space of intermediate states. To this end, we have 
kept fixed the set of SP energies, using the one obtained from the effective 
hamiltonians with the largest number of intermediate states, and calculated 
the energies of the yrast states in 6 Li relative to 4 He, using effective two-body 
interactions V eS that correspond to different values of N max . 

The results are shown in Figs. H] and for the N 3 LO and N 3 LOW poten- 
tials respectively. 




Figure 4: Theoretical energies of 6 Li yrast states relative to 4 He, obtained with the N 3 LO 
potential, as a function of iV max (see text for details). 



Comparing these figures with the previous Figs. [2] and |3] it can be seen 
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2 4 6 8 10 12 14 16 18 

Figure 5: Same as in Fig. 0]for N 3 LOW. 

that a faster convergence is obtained. This may support the choice of smaller 
intermediate state spaces when the SP energies are not derived theoretically 
but taken, for instance, from experiment. 

From now on, our calculations will refer to an intermediate state space 
with N max = 22 and N max = 18 for the N 3 LO and N 3 LOW potentials, 
respectively. 

It is worth noting that the use of a second-order Q-box does not imply a 
loss of generality of the foregoing discussion, since both at second and third 
order in the Q-box expansion the Q space is spanned by the same particle- 
particle, 3-particle 1-hole, and 4-particle 2-hole excitations. Therefore the 
conclusions drawn previously still hold for an effective hamiltonian derived 
including diagrams up to third order in the Q-box expansion. 

Concluding this subsection, we would like to point out that an alternative 
approach to the study of the convergence with respect to the intermediate- 
state space could be to sum at all orders the ladder diagrams of the per- 
turbative expansion by calculating the Brueckner reaction-matrix G of the 
N 3 LO potential. 

However, this is not an easy task, since to manage the operator Q 2p = 
1 — P in the integral equation that defines the G matrix, 
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G(u) = V NN - V NN Q 



2p 



H, 



-Q 2p G(u) 



(40) 



would mean to handle matrix elements of the NN potential Vnn i* 1 a very 
large space spanned by the HO wavef unctions. 

A way to simplify the calculation of the G matrix is to resort to the so- 
called Gt matrix 33( for which plane waves are used as intermediate states: 



Gt(w) — Vnn — VnnQ 



1 



^, n m • ( 41 ) 

^ — Hip 1 V2p 

A main advantage of the Gt matrix is that an exact treatment of the 
projection operator Q 2p can be done employing the Tsai-Kuo method [34J . 

By renormalizing the N 3 LO potential with the Gt matrix procedure, we 
have calculated the energies of the yrast states in 6 Li relative to 4 He, using 
effective hamiltonians derived at the second order in Hi, as a function of 
N m ax- The energies are reported in Fig. [6] (case (a)), and compared with the 
results previously shown in Fig. |2] (case (b)). 




Figure 6: Theoretical energies of 6 Li yrast states relative to 4 He, obtained with the Gt 
matrix derived from the N 3 LO potential, as a function of iV max (see text for details). 
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From inspection of Fig. EJ it is evident that the renormalization of the 
N 3 LO potential strongly reduce the dependence on the number of interme- 
diate states of the perturbative expansion of H c s, and consequently a faster 
convergence with respect to N max is reached. One has to keep in mind, 
however, that, since the Brueckner reaction matrix is energy dependent, the 
results (a) are slightly dependent on the choice of the starting energy e in 
Eq. (USD. 

3.3. Order-by- order convergence 

Now, it is time to focus our attention on the dependence of the effective 
hamiltonian on the order at which the perturbative expansion of the Q-box 
is arrested. 



Historically, this problem was faced first in Ref. [35|, where some selected 
third-order and few fourth-order terms in the G matrix were calculated, and 
the convergence order-by-order of the perturbation series was investigated. 

As mentioned in Subsec. 13.14 we have derived effective hamiltonians 
using Q-boxes at second order (i?|f d ) and third order (H^ d ) in perturbation 
theory. Besides, we have also derived effective hamiltonians calculating the 
Pade approximant [2|1] jii], 26[ of the Q-box (ifpf dfc ). The matrix elements 



of the latter can be found in the Tables in Appendix A 
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Figure 7: Theoretical energies of 6 Li yrast states relative to 4 He, obtained with iff^, H^ dl 
H^f d , and Hpf db derived from the N 3 LO potential (see text for details). In the right side 
of the figure, where an expanded scale is adopted, the Hff t results are omitted. 
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Figure 8: Same as in Fig. [7] for N 3 LOW. 



Actually, in the definiton of the Pade approximant [2|1] (see Eq. f j39l) ) 
also the contribution at first order in Hi comes into play, so it is worth to 
consider also the results obtained using only the diagram of Fig. [1] for the 
one-body hamiltonian, plus the NN bare potential (H°f t ). 

In Figs. [7] and [HI we report the energies of 6 Li yrast states with respect 
to 4 He, obtained with the N 3 LO and N 3 LOW potentials, respectively, and 
calculated with Hfg, Ff d , # 3 f d , and #f ade - 

The large difference between the results with Hff t and those obtained 
taking into account the correlations with the core and the shells above the 
model space (H^ d , H^ d ) evidences the poor description provided by the bare 
NN potential without any renormalization due to long-range correlations. 

Comparing the results shown in the right sides of both figures, we ob- 
serve that the differences in the energies obtained with H^ d and # 3 f d arc 
less notable using the low-momentum N 3 LOW potential. This traces back 
to the fact that the N 3 LO potential, even though its gaussian cutoff is 
around 2.5 fm~ , exhibits some repulsion in the intermediate range, while 
the N 3 LOW potential is a real low-momentum potential. 

This is testified by the fact that in nuclear matter the N 3 LO potentials 
saturates, while the N 3 LOW potential does not; in Fig. [9] we report the 
energy per nucleon in symmetric nuclear matter obtained from a Brueckner- 
Hartree-Fock calculation. Obviously, the larger repulsive components of the 
N 3 LO potential influence negatively the order-by-order convergence. 
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Figure 9: Energy per nucleon in symmetric nuclear matter obtained from a Brueckncr- 
Hartree-Fock calculation with N 3 LO (smooth curve) and N 3 LOW (dot-dashed curve). 



It is worth pointing out, however, that both potentials lead to results 
obtained with H^ db that are very close to those with H^ d , thus supporting 
the hypothesis of a weak dependence of the results on higher-order Q-box 
perturbative terms. 

On the above grounds, in the following discussions we will employ only 
effective hamiltonians derived by calculating the Pade approximant [2|1] of 
the Q-box. 

Finally, it is appropriate to draw attention to the papers by Hjorth- Jensen 



et al. [36l . I371 l38l . |39| where, within the framework of the folded-diagram 
theory, the inclusion of the third-order diagrams in the G-matrix in the cal- 
culation of the Q-box was studied and the order-by-order conver gen ce of 



39 



was 



the effective interaction was examined. The main finding of Ref. 
that the effects of third-order contributions in the T = 1 channel are almost 
negligible. 

3.4- Dependence on the harmonic-oscillator parameter 

As we have mentioned before, the shell-model effective hamiltonian is 
derived using the harmonic-oscillator potential as the auxiliary potential U. 

Since we truncate the number of intermediate states and arrest the per- 
turbative expansion of the Q-box at a certain order, this introduces a depen- 
dence on the value of the harmonic-oscillator parameter. 
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Figure 10: Theoretical energies of 6 Li yrast states relative to 4 He, obtained with the N 3 LO 
potential, as a function of huj (see text for details). 

In order to study this dependence, we have derived from the N 3 LO poten- 
tial three effective hamiltonians using Hcu— 18, 19, and 20 MeV, respectively. 

In Fig. [TUJ the theoretical energies of the yrast states in 6 Li are reported 
as a function of Hco; case (a) refers to effective hamiltonians derived including 
in the Q-boxes diagrams up to second order in perturbation theory. Case (b) 
refers to effective hamiltonians derived including all third-order diagrams in 
the Q-box, and then calculating its Pade approximant [2|1]. 

From inspection of Fig. [TO], it can be observed that while the eigenvalues 
of the second-order effective hamiltonians retain a significative dependence 
on the harmonic-oscillator parameter, the effective hamiltonians derived cal- 
culating the Pade approximant [2|1] of the Q-box are far less dependent on 
fkd. It is worth to note, however, that in case (a) the relative spectra of 6 Li 
are almost independent of ftw. 

The results shown in Fig. [10] highlight the need to include higher-order 
terms in the perturbative expansion of the Q-box. In this connection, for 
the sake of completeness we compare the results obtained in case (b) with 
those obtained without including second- and third-order (V — f/)-insertion 
diagrams in the Q-boxes, i.e. taking into account only the first-order one 
(see Fig. [T|). In Fig. [HI the theoretical energies of the yrast states in 6 Li 
obtained with this procedure (case (c)) are reported and compared with (b). 
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Figure 11: Same as in Fig. [TU](see text for details). 



The results reported in Fig. [TT] show that without including consistently 
order by order the (V — LQ-insertion diagrams in the Q-box, a relevant de- 
pendence on foul may be introduced. 



Finally, it should be mentioned that in the works in Refs. |40|, |4l|], the 
role of (V — C/)-insertion diagrams has been investigated. 

3.5. An application of the Z-box graphical method: (0 + 2) fiw shell-model 
calculations 

Starting from a HO unperturbed hamiltonian, it is natural to extend our 
calculations for the p-shell model space to a (0 + 2) fko one, which includes 
the s<i-shell orbitals too. In such the unperturbed hamiltonian yields 

a non-degenerate model space, the p and sd orbitals being separated by 
AE = hw. 

In Subsec. 12.31 it has been shown that the Z-box graphical method pro- 
vides a simple way to derive effective hamiltonians for non-degenerate model 
spaces. We have therefore found it interesting to employ it to obtain a real- 
istic effective shell-model hamiltonian for the psd model space. 
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To this end, we have calculated the Pade approximant [2|1] of the Q-box 
as a function of the energy e, starting from the chiral N 3 LO potential. The 
Q(e) has been then employed to calculate the vertex-function Z(e), as defined 
by Eq. ( 1321) . The Z(e) is the building block of the graphical method we have 
described in Subsec. 12.31 to derive H cG , whose equations have been solved 
using the secant method algorithm (see Fig. fl~2]) . 

We have to recall now that, when dealing with a (0 + 2) hui model space, 
the results of the diagonalization of the shell-model hamiltonian are affected 
by the spurious center-of-mass motion [42]. The procedure we have followed 
to separate in energy the excitations due to the internal degrees of freedom 
from those with spurious center-of-mass components is the one suggested by 



Gloeckner and Lawson in Ref. 43 



24 



CD 




(a) 0.5 1.0 1.5 2.0 2.5 Offco 



Figure 13: Theoretical energies of 6 Li yrast states, obtained with the N 3 LO potential, 
relative to 4 He as a function of parameter a = log (3. (a) are results with (3 = (see text 
for details). 

Thus we have diagonalized the modified shell-model hamiltonian H' 

H' = H cS + Hp , (42) 

where Hp is /3 times the center-of-mass excitation energy of the A-nucleon 
system 

The spurious components are pushed up in energy by increasing the pa- 
rameter /3, so that one can assume that the low-energy spectrum is free from 
the above components. Usually, the parameter (3 is expressed as a power of 
ten (3 = 10°, and the hamiltonian H' is diagonalized increasing a until the 
low-energy eigenvalues are stable. 

In Fig. [131 we report our results for the yrast states of 6 Li as a function 
of the a parameter, and compare the positive parity spectrum with the one 
reported in Fig. [7| corresponding to the Hcu p-shell model space, using the 
LS method where the Pade approximant of the Q-box has been calculated. 
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From inspection of Fig. [T21 it can be observed that the results are quite 
stable for a > 1.5. Moreover, it should be noticed that the yrast 4 + , 5 + , 
1~, and 4~ states have not been reported since they turn out to be strongly 
affected by center-of-mass spuriosity. 

Finally, it is worth to note that yrast + , 1 + , 2 + , and 3 + are in an excellent 
agreement with those calculated within the ftw model space. 



4. Comparison of realistic shell model with ab initio calculations 

In the previous subsections we have shown how the approximations in- 
volved in the derivation of realistic shell-model effective hamiltonians may 
be kept under control by way of some convergence checks. 

However, to study the accuracy of these approximations one should com- 
pare the final results with those provided by an approach that gives an "ex- 
act" solution of the Schrodinger equation ([I]). 

To this end, we will compare in this subsection our results obtained for 
the p-shell nuclei starting from the N 3 LO potential with those provided from 



the ab initio no-core shell model (NCSM) [15|, [10 . 

All the results shown in the previous subsections have been obtained 
starting from a A-body hamiltonian which is not translationally invariant: 

^E^+EC ■ («) 

i=l i<j=l 

In order to compare realistic shell model with NCSM we have to employ a 
purely intrinsic Hamiltonian, so we have to remove the center of mass kinetic 
energy from the hamiltonian of Eq. 
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Figure 14: Theoretical energies of 6 Li yrast states relative to 4 He, obtained with N 3 LO 
potential, (a) shell-model calculation with an effective hamiltonian derived from Eq. (|44|) . 
(b) shell-model calculation with an effective hamiltonian derived from Eq. (I45|) . (c) NCSM 
calculation. 



We rewrite the intrinsic hamiltonian as follows: 



Hi 



iiit 



(l a) Ei=l 2m + EiO=l (Vi 



ij 



E 



PL 

i=l V 2m 



+ 



4 (\/NN 
l\ v ij 



i<3 



Ui 



_ Pi-Pj \ _ 

Pi _ Pi Pj 

2mA mA < 



(45) 



In Fig. the calculated energies of the yrast states in 6 Li relative to 
4 He are reported. Results labelled with (a) refer to a shell-model calculation 
with an effective hamiltonian derived from Eq. (1441) . while the spectrum 
(b) corresponds to an effective hamiltonian derived from the translationally 
invariant hamiltonian of Eq. (f45|) . The NCSM spectrum is obtained using 
the value of the binding energy as calculated in Ref . [lo[ , combined with the 



calculated excitation energies reported in Ref. [15], and is calculated with 
respect to the 4 He ground state energy provided by the N 3 LO potential 44 . 

From inspection of Fig. [TH it is evident the need to employ a purely 
intrinsic hamiltonian when dealing with light nuclei. Moreover, it can be 
observed that the agreement between the shell-model results and the NCSM 
ones is quite good. 
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Figure 15: Ground-state energies for N = Z nuclei with mass 6 < A < 12. 

At this point it is worth to remind that our shell-model effective hamil- 
tonian is derived for a system with only two valence-nucleons. This implies 
that using this hamiltonian to study nuclei with n > 2 valence nucleons, one 
introduces an approximation due to neglect the 3-, 4-, .. n-body components 
that arise in Hf, even if the original hamiltonian contains only a two-body 
force. 

With this remark in mind, we compare in Fig. [15] the ground-state ener- 
gies, relative to 4 He, for the N = Z nuclei with mass 6 < A < 12 calculated 
within the framework of the realistic shell model (dot-dashed line) and of 
the NCSM (dotted line), with the experimental ones (continuous line) jijj ]. 
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The shell-model ground state energy for 6 Li is close to the one derived with 
the NCSM as already shown in Fig. [HJ while the discrepancy between the 
calculated values becomes more significant going from 10 B to 12 C. This can 
be related to the many-body (> 2) components of Hf s , whose role grows 
with the number of valence nucleons. 

Fig. [in] shows the theoretical RSM and NCSM relative spectra of 10 B 
compared with the experimental one [46|. From this figure it can be seen 
that the RSM excitation energies are in good agreement with the NCSM 
ones, the largest discrepancy being less than 1 MeV. 




RSM NCSM Expt 



Figure 16: Theoretical and experimental spectra for 10 B. The theoretical energies have 
been obtained using the N 3 LO potential within a realistic shell-model framework (RSM) 
and via a NCSM calculation (NCSM). 



It should be noted that in Ref. [lOj it has been pointed out how the 
contribution of a three-body force improves significantly the reproduction by 
the theory of both ground state and excitation energies. 

5. Summary and remarks 

This paper has been devoted to present the current status of the deriva- 
tion of a shell-model effective hamiltonian within the perturbative approach, 
starting from a realistic NN potential. 
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As we have already emphasized, the core of this approach is the perturba- 
tive expansion of the well-known Q-box vertex function. Therefore, we have 
focused our attention on the problems arising from the calculation of the 
Q-box perturbative series. We have tackled different questions, all of them 
dealing with the convergence properties (truncation of the intermediate-state 
space, order-by-order convergence, dependence on the HO parameter). 

At present, it is practically unfeasible to calculate contributions to the 
perturbative expansion beyond the third order ones. Thus we have presented 
results obtained by using the theory of the Pade approximants, starting from 
two different NN chiral potentials that are characterized by different offshell 
properties. Comparing the results obtained with effective hamiltonians de- 
rived calculating the Q-box at third order in perturbation theory and those 
obtained calculating the corresponding [2|1] Pade approximant, we conclude 
that modern chiral NN potentials provide effective hamiltonians scarcely 
dependent on higher-order contributions. 

The foregoing conclusion is also borne out by the comparison of shell- 
model results obtained with an effective hamiltonian derived from the chiral 
N 3 LO potential with those provided by ab initio NCSM calculations for some 
p-shell nuclei. 

Finally, we have applied the recently introduced Z-box graphical method 
to derive a realistic effective hamiltonian for the (0 + 2) fvjj psd model space. 
The results show the versatility of this new method, which may be easily 
applied to non-degenerate model spaces. 

We hope that the discussions and results given in this paper may provide 
useful guide to calculating shell-model realistic effective hamiltonians. 
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Appendix A. Effective hamiltonians 



Table A.l: Single-particle energies (in MeV) of different effective hamiltonians derived 
from the N 3 LO potential and calculating the Pade approximant [2|1] of the Q-box. Labels 
LS, Z-box, and A = 6 indicate SP energies obtained using the LS iterative technique, the 
Z-box graphical method, and starting from a purely intrinsic hamiltonian for A = 6 (see 
Sec. HJ), respectively. 



nlj 


T z 


LS 


Z-box 


A = 6 


0P3/2 


1/2 


4.593 


4.543 


2.747 


0pi/2 


1/2 


7.449 


7.512 


5.886 


0P3/2 


-1/2 


3.915 


3.904 


2.052 


OPl/2 


-1/2 


7.077 


7.163 


5.720 



Table A.2: Same as in Table EH for the N 3 LOW potential. 



nlj 


T z 
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Z-box 


A = 6 


0P3/2 
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4.941 
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3.083 


0pi/2 
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7.647 
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5.655 


0P3/2 


-1/2 


4.001 


3.987 


2.126 


0pi/2 


-1/2 


6.954 


6.892 


5.186 
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Table A. 3: Two-body matrix elements (in MeV) of different effective hamiltonians derived 
calculating the Pade approximant [2|1] of the Q-box from the N 3 LO potential. They are 
antisymmetrized, and normalized by a factor 1/^/(1 + Sj a j b )(l + 3j c j d )- Labels are the 
same as in Table POI 



n,„l„i„ nulhi^ n„l„i, njljij 
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-3.804 


u i'l /2 u i^l/2 Vfl/2 u l J l/2 


o 


1 


0.727 


0.725 


0.358 


u l J 3/2 vfl/2 u l J 3/2 


1 


1 


0.611 


0.615 


0.425 


u f3/2 ".F3/2 ".F3/2 ".F3/2 


2 


1 


-0.879 


-0.876 


-1.056 


u l J 3/2 u f3/2 u f3/2 ^fl/2 


2 


1 


-1.677 


-1.704 


-1.737 


Dt)q /o Ofh /o Dt)q /o Ofh /o 
KJ l J 3/2 KJ l J l/2 w / y 3/2 w / y l/2 


2 


1 


-1.824 


-1.887 


-2.163 


".F3/2 ".F3/2 u f3/2 u f3/2 


o 


-1 


-2.746 


-2.742 


-3.062 


".F3/2 "/ y 3/2 ".Fl/2 ".Fl/2 


o 


-1 


-3.754 


-3.760 


-3.958 


vyl/2 vpi/2 vjj\/2 yjy\l2 


n 


_1 


021 


-0 030 

\J .WOW 


-0 400 

W . t: W W 


"-AF3/2 1J / y l/2 "^3/2 u.Fl/2 


1 


_1 


135 


0.111 


-0 054 

W . W ' ./ t: 


KJ f3/2 KJ f3/2 u f3/2 u i J 3/2 


2 


-1 


-1.367 


-1.370 


-1.564 


u f3/2 u f3/2 u l J 3/2 u l J l/2 


2 


-1 


-1.747 


-1.783 


-1.832 


u l J 3/2 vfl/2 u l J 3/2 ^fl/2 


2 


-1 


-2.205 


-2.295 


-2.538 


u f3/2 ".F3/2 ".F3/2 ".F3/2 


o 


o 


-2.656 


-2.642 


-2.976 


u l J 3/2 u f3/2 u l J l/2 ^fl/2 


o 


o 


-3.751 


-3.809 


-3.932 


KJ f\j2 KJ f\j2 ^yi/2 w / y l/2 


o 


o 


0.121 


0.085 


-0.271 


VfS/2 ' J J J 3/2 ' J J J 3/2 u f3/2 


1 


o 


-1.902 


-1.889 


-2.366 


vy3/2 vps/2 ".F3/2 1J / y l/2 


1 


n 


3 615 


3 659 

'J . w w u 


3 680 

• ) .WOW 


u l J 3/2 w / y 3/2 "/-'1/2 ' J l J 3/2 


1 


o 


-3.581 


-3.619 


-3.698 


KJ I J 3/2 KJ l J 3/2 KJ l J l/2 KJ l J l/2 


1 


o 


2.794 


2.832 


3.072 


0^3/2 Opi/2 0p 3/2 QPj/2 


1 





-3.046 


-3.085 


-3.352 


0^3/2 0j?i/2 Opi/2 0p 3 / 2 


1 





3.224 


3.301 


3.347 


0^3/2 Opi/2 0j?i/2 Opi/2 


1 





1.096 


1.128 


1.153 


Opi/2 0^3/2 Qpi/2 0p 3/2 


1 





-3.109 


-3.209 


-3.422 


Opi/2 0^3/2 Qpi/2 Qpi/2 


1 





-1.143 


-1.205 


-1.247 


Opi/2 Opi/2 Qpi/2 Qpi/2 


1 





-2.113 


-2.166 


-2.850 


0^3/2 0p 3/2 0p 3/2 0p 3/2 


2 





-1.369 


-1.370 


-1.554 


0^3/2 0p 3/2 0p 3/2 Opi/2 


2 





-1.230 


-1.252 


-1.286 


0^3/2 0p 3/2 Opi/2 0p 3/2 


2 





1.207 


1.227 


1.285 


0^3/2 Opi/2 0p 3/2 Opi/2 


2 





-4.280 


-4.331 


-4.760 


0^3/2 Opi/2 Opi/2 0p 3/2 


2 





-2.062 


-2.134 


-2.281 


OPI/2 0p 3/2 Opi/2 0p 3/2 


2 





-4.373 


-4.423 


-4.876 


0^3/2 0p 3/2 0p 3/2 0p 3/2 


%4 





-5.137 


-5.118 


-5.565 



Table A.4: Same as in Table [O] for the N 3 LOW potential. 



Jaja n b l b j b n c l c j c n d l d j d 


J 


T 


LS 


Z-box 


A = 6 


0P3/2 0p 3 /2 0p 3/2 0p 3/2 





1 


-2.828 


-2.815 


-3.221 


0^3/2 0P3/2 OPl/2 OPl/2 





1 


-3.554 


-3.585 


-3.800 


OPl/2 OPl/2 OPl/2 OPl/2 





1 


-0.425 


-0.394 


-0.894 


0^3/2 OPl/2 0p 3 / 2 Opi/2 


1 


1 


0.478 


0.486 


0.291 


0p 3 /2 0p 3/2 0p 3/2 0p 3 /2 


2 


1 


-0.896 


-0.899 


-1.111 


0P3/2 0p 3/2 0p 3 / 2 0j?i/2 


2 


1 


-1.893 


-1.906 


-1.971 


0p 3 /2 Opi/2 0p 3/2 0j?i/2 


2 


1 


-2.077 


-2.109 


-2.447 


0p 3 /2 0p 3/2 0p 3/2 0p 3/2 





-1 


-3.573 


-3.554 


-3.892 


0P3/2 0P3/2 Qpi/2 Opi/2 





-1 


-3.779 


-3.788 


-3.995 


Opi/2 Opi/2 Qpi/a Opi/a 





-1 


-1.033 


-1.014 


-1.477 


0p 3 /2 Opi/2 0p 3/2 Opi/2 


1 


-1 


0.042 


0.042 


-0.161 


0p 3 /2 0p 3/2 0p 3 / 2 0p 3 /2 


2 


-1 


-1.399 


-1.400 


-1.620 


0p 3 /2 0p 3/2 0p 3/2 Opi/2 


2 


-1 


-1.985 


-1.987 


-2.064 


0p 3 /2 Opi/2 0p 3/2 QPj/2 


2 


-1 


-2.478 


-2.498 


-2.791 


0p 3 /2 0p 3/2 0p 3/2 0p 3/2 








-3.523 


-3.508 


-3.862 


0p 3 /2 0p 3/2 QPj/2 Opi/2 








-3.743 


-3.759 


-3.957 


Opi/2 Opi/2 0j?i/2 0j?i/2 








-0.987 


-0.966 


-1.436 


0j5 3 /2 0p 3/2 0p 3/2 0p 3/2 


1 





-2.283 


-2.254 


-2.676 


0p 3 /2 0p 3/2 0p 3/2 Opi/2 


1 





3.631 


3.626 


3.682 


0p 3 /2 0p 3/2 Qpi/2 0p 3/2 


1 





-3.542 


-3.536 


-3.661 


0P3/2 0P3/2 Qpi/a Qpi/2 


1 





2.737 


2.759 


3.130 


0p 3 /2 Opi/2 0p 3/2 Opi/2 


1 





-3.462 


-3.439 


-3.709 


0p 3 /2 Opi/2 Opi/2 0p 3/2 


1 





3.486 


3.464 


3.520 


0p 3 /2 Opi/2 PPj/jj Opi/2 


1 





1.183 


1.197 


1.098 


Opi/2 0p 3 /2 Opi/2 0p 3/2 


1 





-3.444 


-3.423 


-3.656 


Opi/2 0p 3 /2 Opi/2 Opi/2 


1 





-1.273 


-1.288 


-1.254 


0^1/2 Opi/2 0j?i/2 Opi/2 


1 





-3.277 


-3.264 


-3.918 


0p 3 /2 0p 3/2 0p 3/2 0p 3 / 2 


2 





-1.396 


-1.400 


-1.616 


0p 3 /2 0p 3/2 0p 3/2 Opi/2 


2 





-1.393 


-1.396 


-1.453 


0p 3 / 2 0p 3/2 Qpi/a 0p 3/2 


2 





1.372 


1.375 


1.461 


0p 3 /2 OPl/2 0P3/2 Qpi/a 


2 





-4.678 


-4.693 


-5.098 


0p 3 /2 Opi/2 0j9i/2 0p 3/2 


2 





-2.055 


-2.041 


-2.285 


Opi/2 0p 3 /2 Opi/2 0p 3/2 


2 





-4.675 


-4.690 


-5.090 


0p 3 /2 0p 3/2 0p 3/2 0p 3/2 


3 





-5.569 


-5.568 


-6.014 
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Appendix B. Diagrammatics 




Figure B.17: One-body second-order diagrams. The asterisk indicates non-symmetric 
diagrams, which occur always in pairs giving equal contributions. For the sake of simplicity 
we report only one of them. 




Figure B.18: Two-body diagrams up to second order in perturbation theory. For the sake 
of simplicity, for each topology we report only one of the diagrams which correspond to 
the exchange of the external pairs of lines. 
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Figure B.19: Same as Fig. IB. 171 but for third-order diagrams. 
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